SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

COMPOSITES 

STRUCTURAL PROPERTIES 

SURFACES, INTERFACES AND 
THIN FILMS 

ATOMISTIC MODELS 

Received 
14 February 2014 

Accepted 
7 March 2014 

Published 
27 March 2014 



Correspondence and 
requests for materials 
should be addressed to 
X.-Y.L. (xyliu@lanl.gov) 



Revisiting the AI/AI2O3 Interface: 
Coherent Interfaces and Misfit 
Accommodation 

Ghanshyam Pilania 1 , Barend J. Thijsse 2 , Richard G. Hoagland 1 , Ivan Lazic 2 , Steven M. Valone' 
& Xiang-Yang Liu 1 

'Materials Science and Technology Devision, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA and, 
2 Department of Materials Science and Engineering, Delft University of Technology, Mekelweg 2, 2628 Delft, The Netherlands. 



We study the coherent and semi-coherent Al/or-Al 2 0 3 interfaces using molecular dynamics simulations with 
a mixed, metallic- ionic atomistic model. For the coherent interfaces, both Al-terminated and O-terminated 
nonstoichiometric interfaces have been studied and their relative stability has been established. To 
understand the misfit accommodation at the semi-coherent interface, a 1 -dimensional (ID) misfit 
dislocation model and a 2-dimensional (2D) dislocation network model have been studied. For the latter 
case, our analysis reveals an interface dislocation structure with a network of three sets of parallel 
dislocations, each with pure-edge character, giving rise to a pattern of coherent and stacking-fault-like 
regions at the interface. Structural relaxation at elevated temperatures leads to a further change of the 
dislocation pattern, which can be understood in terms of a competition between the stacking fault energy 
and the dislocation interaction energy at the interface. Our results are expected to serve as an input for the 
subsequent dislocation dynamics models to understand and predict the macroscopic mechanical behavior of 
Al/or-Al 2 0 3 composite heterostructures. 



Aluminum is one of the world's most widely used metals, in large part due to its superior strength-to-weight 
ratio, while alumina ( A1 2 0 3 ) is one of the most strongly bonded compounds in existence (with enthalpy of 
formation of 1674.4 kj/mol) 1,2 and is widely used as a catalyst support, for structural ceramics, and as a 
substrate for film growth 3 6 Among the various polymorphs of alumina, corundum (a-Al 2 0 3 ) is the most stable 
phase at the ambient conditions and has been extensively studied 7 The Al/a-Al 2 0 3 interface is of enormous 
scientific and technological significance due to the crucial role it plays in a range of important applications such as 
metal-ceramic composites, protective coating for Al, casting and smelting processes, microelectronics, corrosion/ 
wear protection and catalysis 811 As a consequence, large amount of experimental 8,12 " 22 and theoretical 23 * 33 
research work has been devoted to this interface. 

Based on high-resolution transmission electron microscopy, the primary orientational relationship at the Al/a- 
A1 2 0 3 interface has been observed to be one that matches the close-packed planes and directions in the two 
phases 17 In this orientation the Al(lll) plane is parallel to the Al 2 O 3 (0001) basal plane and the following 
directions are parallel to each other: 



(101)^11(0110), 
(121)^11(2110) 



A1 2 0 3 



(1) 

(2) 



As is the case for other similar metal/ A1 2 0 3 interfaces 2,34,35 , the stoichiometry of clean Al/a-Al 2 0 3 interface is a 
function of the oxygen partial pressure. Previous first principles based analysis 25 has suggested that the clean Al/ 
A1 2 0 3 interface can be stoichiometric, oxygen terminated, or aluminum terminated, depending on the ambient 
oxygen partial pressure. For the relatively low oxygen partial pressures 36 , usually associated with sessile drop 
experiments, the interface is predicted to be stoichiometric. However, for relatively high oxygen partial pressures, 
prevalent during fracture experiments, the interface has been suggested to be oxygen terminated 25 . Recent 
experimental 37 and theoretical 31 studies reporting atomically abrupt Al(liquid)/a-Al 2 0 3 (solid) interface stabilized 
by self-regulated interfacial Al vacancies have suggested an aluminum terminated interface. 
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Electronic structure and work of adhesion have been obtained 
from first principles based quantum mechanical computations by 
Batyrev and Kleinman 26 and Siegel et al. 2S,2 ' > assuming a coherent 
interface. Streitz and Mintmire 32,33,38 employed a force field account- 
ing for the electrostatic interaction associated with charge transfer in 
their MD simulations to study the solid A1/A1 2 0 3 interface. However, 
unlike the experimentally-observed abrupt interface 17 , the force field 
predicted that O atoms should rapidly diffuse into the Al lattice, 
resulting in a highly disordered region at the interface 33 . Zhang 
et a\? a have studied the adhesion and nonwetting to wetting trans- 
ition in the Al/a-Al 2 0 3 interface using a reactive force field (ReaxFF) 
approach 3 '. They find that the evaporation of Al atoms and diffusion 
of O atoms in a- A1 2 0 3 primarily lead to the wetting of liquid Al on 
the oxide surface. 

Nanolayered laminated composites of aluminum and aluminum 
oxide with interlayer spacing ranging from 50 to 500 nm have been 
prepared in the past 40 and found to exhibit interesting deformation 
mechanism with extensive ductility. However, an atomistic under- 
standing of the observed mechanical behavior is presently lacking. 
Although available dislocation dynamics models are capable of 
exploring evolution of interface dislocation networks and can predict 
the microscopic properties of the nanolayered composites during 
mechanical deformations at large length scales, such models depend 
upon atomistic input to define interfacial properties and dislocation 
reaction rules. A deeper understanding of the interface structure can 
be crucial for understanding the macroscopic mechanical behavior of 
the nanolayered composite heterostructures. 

However, despite the significance of interfacial structure, studies 
addressing them have remained quite scarce in the literature. The 
reason for this scarcity is two fold. On the one hand, accurate ab initio 
simulation techniques are not efficient enough to practically deal 
with such large systems with several thousand systems 41 . On the 
other hand, sophisticated atomistic potentials required for molecular 
dynamics (MD) simulations, that can treat both metal and ceramic 
systems on equal footing and adequately describe charge transfer 
effects at the interface, were unavailable until recently. As a first step 
towards this direction, the goal of the present study is to elucidate the 
coherent and semi-coherent interface structures at the Al(lll)/o!- 
Al 2 O 3 (0001) interface. This work is expected to provide novel 
information on dislocation patterns that can form at the interface 
to accommodate misfit strain as well as on the factors that govern the 
formation of these patterns. The resulting knowledge from the pre- 
sent study is expected to serve as an input for the subsequent dis- 
location dynamics models 42-43 to understand and predict the 
macroscopic mechanical behavior of Al/a-Al 2 0 3 composite 
heterostructures. 

Results 

Coherent Interfaces. Employing the aforementioned experimentally 
observed orientational relationship 17 , the Al/a-Al 2 0 3 interface was 
modeled in a slab geometry periodically repeating in the plane of 
the interface. Further details of the aomistic model used and 
the computational methodology employed can be found in the 
Methods section presented towards the end of the manuscript. To 
start with, both the O-terminated and the Al-terminated interfaces 
were considered in our atomistic simulations. Since our computed 
lattice constant for the fcc-Al (a A] = 4.05 A) and the x-alumina 
( fl Al 2 o 3 =4.762 A and cm 2 o 3 — 12.99 A) are in excellent agreement 
with the corresponding experimentally measured 17 values, we 
recover the correct magnitude of the measured misfit strain 



and v^Qa^Oj in the aluminum and alumina crystals, respectively. On 



(5=2 



3/2a A l — flAl 2 o 3 



= 0.041. 



y3/2fl A i-(-fl A ] 2 03 
The crystal geometry in the 2D interface is as follows. The periodic 
lengths in the parallel directions [ 1 0 1 ] A1 and [Ollo]^ q 3 are aai 



the other hand, the periodic lengths in the directions [121] A1 and 
[2H0] Al2 o 3 



are \/-a A i and a A i o 3 in the aluminum and alumina 



crystals, respectively 
aluminum match one 



For a coherent interface three periods in 
period in alumina along the direction 



[101] A1 1 1 [0110] a12 q 3 , while one period in aluminum matches one 
period in alumina along the normal direction (i.e., 
[T2l] A1 || [2110] a12 q 3 ). As shown in Eq. (3), this produces a 
compressive strain of 0.041 in both the directions on the aluminum. 

Figure 1(a) shows our fully relaxed interfacial geometry between 
the Al (111) and the O-terminated a-Al 2 0 3 (0001) at 0 K. We find 
that the computed interface geometry agrees well with the one 
obtained using first principles density functional theory based calcu- 
lations 27 carried out using generalized gradient approximation of 
Perdew et al. 44 We further note that the atomic arrangement in the 
vicinity of the interface plane is strongly influenced by the a-Al 2 0 3 
structure and leads to a splitting of the Al metal layer closest to the 
interface into three distinct sub-layers (namely, layers Ai, A 2 and A 3 , 
as marked in Fig. 1(a)) after relaxation. The Al atoms in the first two 
sublayers (i.e., the layers Aj and A 2 ) occupy the sites that are dictated 
by a natural continuation of the a- A1 2 0 3 lattice. That is, the Al metal 
atoms in the first two sub-layers are pulled toward the O-terminated 
alumina surface, occupying the sites that Al would occupy in 
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Figure 1 | Coherent interface models. Schematic of Front (left) and side 
(right) views of the relaxed structure of (a) O-terminated and (b) Al- 
terminated coherent Al/ y- A1 2 0 3 interfaces. For clarity Al atoms on the 
metal side of the interface have been distinguished (colored differently) 
than that of the AUO, side. 
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bulk alumina. By relaxing in this way, all the dangling bonds on the 
oxygen ions at the interface are saturated. The Al atoms in the third 
sub-layer (i.e., the layer A 3 ), are not directly bonded to the interfacial 
O atoms and occupy a hollow site in a slightly distorted hexagonal 
lattice formed by the first two Al metal sub-layers. 

The relaxed interface geometry between the Al (111) and the Al- 
terminated a-Al 2 0 3 (0001) interface is presented in Fig. 1(b). In their 
first principles investigation, Zhang et al. 2s have found this surface 
configuration as the most stable after considering a number of inter- 
face configurations of Al atoms in the metal layer adjacent to the 
interface with respect to the alumina surface, reached by effectively 
sliding the metal layer along the alumina surface. It is interesting to 
note here that the Al- terminated interface shown in Fig. 1(b) can be 
formed by removing the third sub-layer A 3 (highlighted in Fig. 1(a)) 
of the Al metal atoms in the O-terminated interface. Therefore, the 
chemical potential of the Al atoms in the layer A 3 can be used to 
evaluate the relative stability of the two interfaces. We find that the 
transfer of the Al atoms in the the layer A 3 of the O-terminated 
interface to a bulk Al metal reservoir is indeed fhermodynamically 
favorable, resulting in the lower-energy Al-terminated interface. 
This indicates the superior stability of the Al-terminated interface 
in a low oxygen chemical potential regime. According to our 
CTIP + RFMEAM potential calculations the Al-terminated interface 
is more stable by 1.98 J/m 2 , while parallel DFT-PBE calculations 
also show a qualitative agreement favoring the Al-termination by 
0.90 J/m 2 over the O-termination. 

Based on interface adhesion energy analysis, previous DFT calcu- 
lations have also predicted the Al-termination for the Al(lll)/ 
a-Al 2 0 3 solid-solid interfaces 28 . We also note that Al(liquid)/ 
a-Al 2 0 3 (solid) interface has also recently been found to be atomist- 
ically sharp with Al termination in both in situ high-resolution trans- 
mission electron microscopy experiments 37 as well as through first- 
principles MD simulations 31 . Therefore, we chose the Al-terminated 
interface to study the structure of misfit dislocations at this interface. 
However, we do not expect the pattern of misfit dislocations formed 
at the O-terminated interface to be drastically different than that of 
the Al-terminated interface, since the two interfaces are related to 
one another through addition or removal of an extra sub-layer of Al 
atoms in the vicinity of the interface plane. Such mass transport 
should be easily achievable by a suitable mechanism such as disloca- 
tion glide or diffusion of Al atoms on the metal side of the interface, 
especially at elevated temperatures. 



coherent. To construct the coherent dichromatic pattern (CDP) ref- 
erence state 45 , we compressed and stretched the Al layers on the Al- 
metal side and on the A1 2 0 3 side, respectively, by equal amounts 
along the direction [lOlJ^H [OlIO]^ 0 . The top and side views of 
the reference state are depicted in Fig. 2(b). After constructing the 
reference state, a correspondence list between nearest neighbor Al 
atoms in the metal layer and the ceramic layer is compiled in this 
configuration. Disregistry vectors Ar are finally computed as 

s8 



Ar=rf- 



•i' 



(4) 



where is the relative position between the f* atom and the j" 1 atom 
that form a pair in the reference and r® is the relative position 
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ID Misfit Dislocations. Since the lattice parameters in the Al and the 

a-Al 2 0 3 phases are not identical, a long range elastic stress field 
accompanies a perfectly coherent interface. The coherency field 
can be relieved through formation of a 2D misfit dislocation 
network at the interface. However, before going to the case of a 2D 
network, we first investigate the case of misfit dislocation arrays only 
along the [101] A1 | [OlIO]^ 0 direction, while assuming a perfectly 
coherent interface along the [l2l] A1 || [2110]^ 0 direction. For this 
ID misfit model, the interface supercell dimensions along the 
directions [l0l] A1 || [0110]^ and [M]^ [2110]^ were 65.96 A 
and 4.76 A, respectively. In this case, to minimize overall strain in the 
[101]^ direction, twenty-three aluminum periods were made to 
match eight alumina periods, instead of the 3 : 1 ratio that 
produces a coherent interface. This reduces the overall strain in 
this direction from 0.041 to 1.78X10 -3 . In the other direction the 
period ratio was kept at 1 : 1, leaving the strain unaltered at 0.041, as 
mentioned above for a perfectly coherent interface. 

To quantify the misfit dislocation structure, we measured disreg- 
istry across the interface plane after relaxation. For the disregistry 
measurement the two Al layers nearest to the relaxed interface plane 
(one on A1 2 0 3 side and one on the Al metal side) were taken into 
account (c.f. Fig. 2(a)). To carryout this measurement a reference 
state must be chosen, in which the two adjoining crystals are 
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Figure 2 | ID misfit dislocation model, (a) Relaxed and (b) coherent 
dichromatic pattern (CDP) reference structures containing two Al layers 
nearest to the interface that were used for the disregistry analysis of the ID 
misfit dislocation model. Al atoms on the ceramic and metal sides of the 
interface have been shown in purple and lavender colors, respectively, (c) 
Disregistry vectors Ar computed as the difference between the positions of 
the selected pair of Al atoms in the relaxed and reference configurations, 
(d) Screw and edge components of the disregistry vectors. 
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between the same pair of atoms at the interface 46,47 . The disregistry 
between the two Al atomic planes on either side of the interface is 
thus a position dependent vector field. Furthermore, since the 
Burgers vectors of the misfit dislocation arrays are expected to lie 
within the interface plane, we restrict our attention to disregistry 
components along the directions [101]^ and [12 1] A1 only. We find 
that the resulting disregistry pattern is mainly dominated by the 
homogeneous deformation used in constructing the strained ref- 
erence configuration. As a further analysis, in Fig. 2(d) we plot the 
components Ar c, and Ar h of the disregistry vectors along directions 
t, and h, parallel and normal to the set of misfit dislocation lines 
(laying along the direction [l2l] A1 ), respectively. We find that the 

screw component along the direction £ remains close to zero 
throughout and exhibits some small deviations owing to a slight 
buckling of the Al metal layer adjacent to the interface plane after 
relaxation. On the other hand, the edge component along the dir- 
ection h simply shows a straight line with a constant slope. The shape 
of this pattern is indicative of the presence of a purely Vernier-type 
misfit dislocation with wide overlapping dislocation cores. 

2D Misfit Dislocation Network. For the 2D misfit dislocation model, 
the dimension of the supercell along the directions [l0l] A J 
[0ll0] Al2Oi and [12 T\ M \\ [2110]^ were 65.96 A and 114.29 A, 
respectively. Here, twenty-three aluminum periods were made to 
match eight alumina periods in the [101]^ direction and twenty- 
three aluminum periods were made to match twenty-four alumina 
periods in the [121]^ direction. The dimensions of the simulation cell 
were chosen so as to minimize the strain required for maintaining the 
periodic boundary conditions along the two orthogonal sets of 
directions in the interface plane while maintaining a minimal 
system size. 

After relaxation, the interface structure was analyzed by mea- 
suring the disregistry across the interface plane, as described above 
for the 1 D misfit dislocation model. Disregistry vectors Ar computed 
as the difference between the positions of the corresponding pair of 
Al atoms (in the adjacent layers above and below the interface plane) 
in the relaxed and reference configurations are shown in Fig. 3(a). 
The edge component Ar-h of the disregistry vectors along the 
[l2lO] A j 2 Q 3 direction is plotted in Fig. 3(b), while the screw compon- 
ent remains zero throughout and is not shown in the figure. Similar 
pattern for the disregistry component was also obtained for the 
[2ll0] Al2O3 and [Il20] Al2Oj directions, indicating that there exist 
three sets of edge dislocations at an angle of 120° from each other. 
A closer look at the Fig. 3(b) reveals steps-like flat periodic regions 
which are the result of the tendency of the interface deformation to 
localize into compact coherent regions, while the steepest slope 
correspond to the dislocation cores. Several parameters character- 
izing the dislocation network can be read from the Fig. 3(b). 
Difference between successive horizontal flat regions of the disreg- 
istry component Ar-h is the magnitude of the net edge component b e 
= 2.48 A of the Burgers vectors of the dislocations. On the other 
hand, the spacing s = 57.1 A between the parallel sets of dislocations 
along the [I2I0] A j 2 q 3 direction is the distance between the successive 
horizontal regions along the abscissa. 

The 2D dislocation network for the relaxed Al/a-Al 2 0 3 interface 
as predicted from our disregistry analysis is presented in Fig. 3(c). 
The interface structure reported is also found to be consistent with 
the Frank-Bilby equation 45,48 50 , details of which are provided in 
the Method section of the present manuscript. For clarity, only three 
Al-layers (i.e., layer A, consisting both the sub-layers A t and A 2 , 
layer B and layer C as depicted in the inset of the Fig. 3(c)) at the 
interface are included in the figure with atoms in each of the three 
layers colored according to their charges. The positions and line 
directions for the identified three sets of edge dislocations are pre- 
sented as solid black lines and black arrows in the figure, while the 



Burgers vectors are depicted as blue arrows. The dislocation network 
clearly reveals coherent (fcc-like) and stacking fault (hep-like) 
regions, separated by the dislocations lines. Some of the coherent 
(hexagonal) and stacking fault (triangular) regions are highlighted 
in green and blue, respectively. Here note that, unlike layers B and C, 
the layer A (i.e., the first Al-layer at the interface on the ceramic side) 
is not a closed pack layer of Al atoms and has a hexagonal lattice 
instead. Therefore, the terms fcc-like and hep-like are used here only 
in a qualitative and local sense with their scope limited to the inter- 
face region alone. 

Figure 4(a) shows the charge distribution of the Al and O atoms 
across the interface along the [lll] Ai || [000l] A j 2Oj direction. It can be 
seen that after a couple of atomic layers away from the interface, on 
both the metal and the ceramic sides, bulk-like behavior in the charge 
distribution is recovered. On the metal side the Al atoms remain 
charge neutral throughout. Surface effects at the free surface on the 
A1 2 0 3 side lead to a deviation in the charges from their respective 
bulk values. The in-plane charge distribution in the first three Al 
metal layers closest to the interface is presented in Figures 4b-d. 
Interestingly, the charges in the first layer B show a charge distri- 
bution pattern in close correlation with the coherent regions of the 
predicted interface structure. This correlation can be attributed to the 
slight variations in the local coordination of Al atoms influenced by 
the interface dislocation structure. As we go deeper in the metal side 
away from the interface, charges quickly diminish in the magnitude 
and the correlated pattern gradually fades. For instance, the layer C 
shown in Fig. 4c has an average charge that is about one order of 
magnitude lower than that of the layer B, nonetheless still shows a 
pattern in the charge distribution. However, the third metal layer 
(i.e., the layer D) further away form the interface (c.f. Fig. 4(d)) 
reveals charges that are already negligible as compared to the other 
two layers and do not carry any correlation with the interface 
structure. 

We note that the interface dislocation structure presented in 
Fig. 3(c) was obtained at 0 K. In the molecular dynamics run, we 
started the simulation at 100 K, and systematically removed the 
kinetic energy of the system in each step during the dynamics. 
Convergence was reached when the temperature was dropped to 
0 K and all the force components on the atoms were lower than 
0.01 eV/A. However, owing to sluggish kinetics at such low tempera- 
tures one may not have converged to the thermodynamically most 
stable state of the system under investigation. To further confirm the 
anticipated structural change in the interface structure at elevated 
temperatures, we carried out annealing of the 0 K relaxed structure at 
500 K for 200 ps, followed by a ramping down to 0 K over a 100-ps 
interval and finally relaxation at 0 K for another 200 ps. 

The final interface structure is presented in Fig. 5(a) with the three 
sets of dislocation lines and their Burgers vectors depicted in the 
figure and is found to be 0.18 J/m 2 lower in energy than the one 
obtained at 0 K. The structure reveals that the annealing treatment 
only has an effect of rigidly shifting the relative positions of any 
of the two parallel sets of dislocations with respect to the third 
set without changing the relative distance within a given set of 
dislocations. This leads to an expansion of the alternate triangular 
hep-like regions (highlighted in blue) around a coherent hexagonal 
region shown in Fig. 3(c), while the other three regions shrink to 
form nodes as can be seen highlighted in green in Fig. 5(a). In other 
words, the hep-like regions transform from equal-area pairs into 
pairs consisting of a large area triangle and a small area triangular 
node. Furthermore, the coherent fcc-like regions transform into 
almost triangular shape regions after the annealing treatment. It 
can also be easily seen that in this interface structure each node is 
surrounded by alternating fcc-like and hep-like regions and has a 
three-fold symmetry in the interface plane. We note that such inter- 
face structures have recently been predicted for interfaces formed 
between various fee metals 51 . 
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Figure 3 | 2D misfit dislocation network, (a) Disregistry vectors in the interface plane computed with respected to a CDP reference structure, (b) Net 
edge component of the disregistry vectors along the [2110]^ 0 direction, (c) The 2D misfit dislocation network for the relaxed AI/2-AI2O3 interface 
predicted at OK. Only three Al-layers (namely layers A, B and C as shown in the inset) are included in the figure with atoms in each of the three layers 
colored according to their charges. 



Discussion 

To gain a deeper insight into the computed dislocation patterns, next 
we investigate various factors that might play a role in controlling the 
predicted interface structure. The energy difference per unit area 
between the fcc-like coherent regions and the hep-like stacking fault 
regions is one of the crucial factors that has a strong influence on the 
observed dislocation pattern. In one extreme limit, when the relative 
energy of the stacking fault regions is much higher than that of the 



coherent regions, a hexagonal Wigner-Seitz mesh-like network of 
misfit dislocations around the points of the best matching (indicated 
as •) in the coherent region is expected 52 , as depicted in the top panel 
of Fig. 6(a). 

The other extreme situation arises when the energy difference 
between the two regions is relatively small, as is the case in the present 
study. We calculate the difference per unit area between the hep-like 
stacking fault regions and the fcc-like coherent regions to be 
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Figure 4 | Charge distributions, (a) The charge distribution profile of the Al and O atoms normal to the interface plane along the [1 1 1] ^ 1 1 [0001]^ 0 
direction for the 2D misfit dislocation network relaxed at OK. The in-plane charge distribution in the (b) first, (c) second and (d) third Al layers closest to 
the interface on the metal side. 



0.03 J/m 2 (compared to our DFT computed value of 0.08 J/m 2 ). In 
this case, the alternate nodes in the hexagonal Wigner-Seitz mesh- 
like network of misfit dislocations (shown as © in the top panel of 
Fig. 6(a)) grow at the expense of the coherent regions until the ratio of 



their interfacial areas approaches unity. As a result, a triangular mesh 
of Shockley partial dislocations, as shown in the bottom panel of 
Fig. 6(a), is observed. Note that going from one extreme case of the 
hexagonal mesh to the other extreme case of the triangular mesh, the 
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Figure 6 | Dislocation network patterns at the A1-A1 2 0 3 interface, (a) 

Schematic illustration depicting the effect of the relative energy difference 
per unit area between the fcc-like coherent regions and the hep-like stacking 
fault regions at the interface. A high stacking fault energy favors a 
hexagonal network of dislocations, while a low stacking fault energy leads 
to a triangular dislocation mesh, (b) and (c) represent the schematics of the 
dislocation patterns obtained in the present study at 0 K and after 
annealing at 500 K, respectively. The dislocation line directions and 
Burgers vectors are shown by black and blue arrows, respectively. See text 
for further details. 

area ratio of the stacking fault regions to the coherent regions varies 
from zero to one. The two interface dislocation patterns, reported in 
Fig. 3(c) and Fig. 5, that we have predicted here represent two inter- 
mediate situations between the two possible extreme cases with the 

ratio of the two regions being - and slightly lower than 1 and are 

depicted schematically in Fig. 6(b) and Fig. 6(c), respectively. 

Dislocation interaction is the second important factor that con- 
trols the observed dislocation pattern at the interface. To further 
elucidate this factor, in Fig. 6(b) we have schematically captured 
the dislocation pattern obtained through structural relaxation near 
the absolute temperature. The dislocation lines and Burgers vector 
directions are shown by black and blue arrows, respectively. The 
center of the coherent region is marked by symbol •, while the sym- 
bols O and © mark the centers of the alternate triangular stacking 
fault regions. It is interesting to note that the dislocation interactions 
are attractive and repulsive in the alternate triangular regions colored 
in green and yellow, respectively, in Fig. 6(b). Therefore, owing to the 
attractive interactions, a high temperature relaxation leads the green 
triangular regions to shrink and collapse in a node. On the other 
hand, the yellow triangles belonging to the second set of the stacking 
fault regions, shown in Fig 6(b), grow larger upon relaxation to push 
the dislocations further apart in order to minimize their repulsive 
interactions. This leads to the dislocation pattern shown in Fig. 6(c) 
(obtained after annealing of the 0 K relaxed structure at 500 K for 
200 ps), while the total dislocation length per unit interfacial area 
does not change. 

In summary, we have studied the coherent and semi-coherent Al/ 
a-Al 2 0 3 interfaces with experimentally known orientation relation- 
ship [110](111) A1 || [l0l0](0001) Al2O3 . Our atomistic simulations 
employed a "reference free" version of the modified embedded atom 
method (RF-MEAM) potential that is capable of adequately describe 



the charge transfer effects at the interface. The simulation results 
predict an interface dislocation structure with a network consisting 
of three sets of parallel dislocations, each with a pure edge character, 
surrounding hexagonal coherent regions and triangular stacking 
fault regions at the interface. The pattern formation can be rationa- 
lized in terms of a competition between the stacking fault energy and 
the dislocation interaction energy and has implications for disloca- 
tion dynamics models aimed at understanding the macroscopic 
mechanical behavior of the Al/a-A^O^ composite heterostructures. 

Methods 

For our atomistic simulations, we employed a recently developed potential for the Al- 
O system 53 , generated following a two-step construction method that is generally 
applicable to ionic systems. As non-electrostatic part of the potential we use a ref- 
erence free version of the modified embedded atom method (RF-MEAM) potential 
that includes angular forces. For the electrostatic part we use the charge transfer ionic 
potential (CTIP) model proposed by Zhou et al. (here after referred to as ZWFN) 54 . 
The modified potential has been shown to achieve essential improvements over the 
ZWFN potential, derived by the incorporation of angular dependent forces and 
angular screening, as well as by the systematic fitting of the potential to a large energy 
database of different Al^O^, crystals, over a range of lattice constants and elastic 
deformations. The current potential yields more realistic atomic charges, crystal 
energies, lattice constants and correct ordering of phase stabilities for a range of Al-0 
systems. Moreover, O-O interactions as well as O interactions with a pure Al (111) 
and Al-terminated a- A1 2 0 3 surfaces are realistically captured within this formulation. 

For the coherent interfaces, aforementioned experimentally observed orientational 
relationship 17 with the Al-terminated interface was employed. The Al/x-Al 2 0 3 
interface was modeled in a slab geometry periodically repeating in the plane of the 
interface. The Al 2 0 3 part of the interface contained six oxygen atomic layers and 
fourteen Al layers, and the metal part contains seven layers. A vacuum region of more 
than 35 A was used between the periodically repeating slabs in the direction normal to 
the interface plane. The final interface structures are determined by relaxing all the 
internal coordinates until all the force components on the atoms were lower than 
0.01 eV/A. For the relaxation of the interface structures we have used an anneal-and- 
quench molecular dynamics (MD) formalism as described in the text, as opposed to 
the molecular statics. 

To compare with the results of atomistic simulations obtained using our 
CTIP + RF-MEAM potential, the density functional theory (DFT) 55,56 computations 
for the coherent interface structures were performed using the Vienna ab initio 
simulation package (VASP) 57 . VASP uses a plane- wave basis set for the expansion of 
the single- particle Kohn-Sham wave functions. The projector- augmented wave 
(PAW) frozen-core potentials 58 ' 59 were used to describe the electron-ion core inter- 
action, and the quantum mechanical part of the electron- electron interactions was 
approximated using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation 
functional 44 . Sampling of the irreducible wedge of the Brillouin zone is performed 
with a regular Monkhorst-Pack grid of special rc-points 60 . To facilitate the numerical 
convergence associated with integrating the 0-K Fermi-Dirac distribution step- 
function, partial occupancies of the single- particle wave functions are introduced 6 ' 62 . 
The electronic wave functions are expanded in plane waves up to a cut-off energy of 
500 eV. A vacuum region of 15 A was used to prevent any spurious unphysical 
interactions between the periodically repeating slabs in the direction normal to the 
interface plane. The final interface structures are determined by relaxing all the 
internal coordinates using the conjugate gradient method 63 until all the Hellmann- 
Feynman forces 64 ' 65 were less than 0.01 V/A. 

The Frank-Bilby equation 45,48 " 50 used in our analysis is given as follows: 

Ebisiny 1 „ , , 

; -^ = Or c p. (5) 

Here, index / runs over all the unique sets of dislocation lines present at the interface, 
each with Burgers vector b\ and separation length l! between the parallel lines within a 
given set. y' denotes the angle between a dislocation line direction if and the probe unit 
vector p, measured in the anti- clockwise direction. D rc represents a deformation 
matrix that maps a unit vector in the coherent reference configuration onto the 
interface structure containing dislocations. 
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